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Abstract 



It has been observed that decay rate for proton emission from ^3/2 single 
particle state is systematically quenched compared with the prediction of a 
one dimensional potential model although the same model successfully ac- 
counts for measured decay rates from si/2 and /iii/2 states. We reconcile this 
discrepancy by solving coupled-channels equations, taking into account cou- 
plings between the proton motion and vibrational excitations of a daughter 
nucleus. We apply the formalism to proton emitting nuclei ^^''■^^^Re to show 
that there is a certain range of parameter set of the excitation energy and 
the dynamical deformation parameter for the quadrupole phonon excitation 
which reproduces simultaneously the experimental decay rates from the 2^3/2, 
3si/2 and I/111/2 states in these nuclei. 
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Typeset using REVT^ 



Physics of nuclei close to the neutron and proton drip lines is one of the most active 
and exciting research areas of the current nuclear physics. Nuclei beyond the proton drip 
line are unstable against proton emission, but, since a proton has to penetrate the Coulomb 
barrier, their lifetime is sufficiently long to study their spectroscopic properties. Thanks to 
the recent experimental developments of production and detection methods, a number of 
ground-state as well as isomeric proton emitters have recently been discovered, which has 
stimulated many experimental and theoretical works P|-pil[]. 

For proton emitters in the A ~ 150 region, proton emissions from the l/iii/2, 3si/2, 
and 2d^/2 orbitals have been observed. It has been pointed out that a spherical calculation 
based upon a one dimensional optical potential with spectroscopic factor estimated in the 
BCS approximation systematically underestimates the measured decay half-lives for proton 
emissions from the 2(^3/2 state, while the same model works well for emissions from the l/iii/2 
and 3si/2 states in both odd-Z even-N nuclei and odd-Z odd-N nuclei [^^. For the ^^^Lu 
nucleus, this discrepancy was attributed to the effects of oblate deformation of the core 
nucleus ^^''Yb, which alter both the decay dynamics and the BCS occupation probability 
of the valence proton P, p!5| , p0| . The coupled-channels calculations with (32 ~ —0.15 have 



successfully explained the measured decay half-lives for this nucleus |T5| . However, for proton 
emitters such as ^^^Ta, ^^°Re, and ^^^Re, the static deformation parameter of the core nuclei 
i55gf^ 159^^^ and ^6°W is estimated to be (32 = -0.053, 0.080, and 0.089, respectively, based 
on the macroscopic- microscopic mass formula |2l|], and thus the deformation effects will 
be much smaller. Those nuclei are nearly spherical, and vibrational excitations should be 



considered instead of the deformation effects and the associated rotational excitations 22 



For vibrational nuclei, the excitation energy of collective excitations is in general signifi- 
cantly larger than that for rotational nuclei. This makes the channel coupling effects weaker 
in spherical nuclei. We notice, however, that the penetration probability is extremely sensi- 
tive to other degrees of freedom especially at deep barrier energies, regardless of the value 
of their excitation energies It is thus of interest and important to explore the role 

of collective core excitations during proton emission decays of spherical nuclei. 

The aim of this paper is to solve the coupled-channels equations for spherical proton 
emitters in order to investigate whether the effects of vibrational excitations of the daughter 
nucleus consistently account for the measured decay half-life from the l/iii/2, 3si/2, and 
2(^3/2 states. We particularly study proton emissions from the l/iii/2 and 3si/2 states of 
^*^^Re and from the 2^3/2 state of ^^°Re 0, assuming that the vibrational properties are 
identical between the core nuclei ^^°W and ^^^W. We shall exclude in our study the proton 
emitter ^^^Ta, since the experimental data is somewhat ambiguous due to the competition 
between the proton emission and the /5'*"/EC decay. Since the properties of excited states 
of these proton-rich nuclei are not known, we shall study the dependence of the decay 
rate on the excitation energy as well as on the dynamical deformation parameter of the 
vibrational phonon excitation of the daughter nucleus. We will then extract a possible 
combination of these two from the experimental data. We will also discuss the dependence 
on the multipolarity of the phonon state. 

We consider the following Hamiltonian for the spherical proton emitting systems: 

H = V + V,oup{r, a) + H„t, (1) 
where r is the coordinate for the relative motion between the valence proton and the daughter 
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nucleus, and /i the reduced mass for this motion, a is the coordinate for the vibrational 
phonon of the daughter nucleus. It is related to the dynamical deformation parameter as 
^^i^ ~ T^^'-'^-^M ~'~ where A is the multipolarity of the vibrational mode and 

a^^(aA^) is the creation (annihilation) operator of the phonon. H^i-^, = huoYl^o^y^^ciXn is the 
Hamiltonian for the vibrational phonon. In this paper, for simplicity, we do not consider the 
possibility of multi-phonon excitations, but include only excitations to the single phonon 
state. The coupling Hamiltonian Vcoupi^.a) consists of three terms, i.e., Vcoup{,f,a) = 
K;L«p(^5 Q;) + Ko«p(^' + K;o«p(^5 The nuclear term reads 

VLZlir, «) = , (2) 

where the dot denotes a scalar product. We have assumed that the nuclear potential is given 
by the Woods-Saxon form. Notice that we do not expand the nuclear potential but include 
the couplings to the all orders with respect to the phonon operator a p5|,p6|. The matrix 



elements of this coupling Hamiltonian are evaluated using a matrix algebra, as in Ref. p5 
As for the Coulomb V}^p as well as the spin-orbit Vj^^^^ terms, the effects of higher order 



couplings are expected to be small [^, and we retain only the linear term. The Coulomb 
term thus reads 

V}^Jp{r, «) = ^ + ^ TTT^ f ^) ■ Ur) (for r > R,) (3) 




^ ax-Yxir) iforr<R,), (4) 



2A + 1 \R, 



C ' 



where is the atomic number of the daughter nucleus and Rc is the charge radius. For 
the spin-orbit interaction, we express it in the so called Thomas form 



(5) 



where /(r) = 1/[1 + exp((r — Rso)/aso)]- The last term in Eq. can be decomposed to a 
sum of angular momentum tensors using a formula 



In order to solve the coupled- channels equations, we expand the total wave function as 
v^,^(r, a) = E E -^"^f^ I iypnl)jm) , (7) 

where 
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{r,a\{lpjpnl)jm) 



E 

nipinf 



{jpmplmi\jm)yj^i^^^ {r)(f)nimi (a) 



(8) 



(j) being the vibrational wave function. We need to compute the couphng matrix elements of 
the operators a\ ■ Tx, where T^^ is either Yx^ or [1a±i (— "i V x a)Y^^^\ These are expressed 



in terms of the Wigner's 6-j symbol as [28 



{ilpjWl')jm\ax ■ Tx\{lpjpnl)jm) 



(9) 



For transitions between the ground and the one phonon states which we consider in this 
paper, the reduced matrix element (0n'/'|l'^A||0n7') is given by /5a- The reduced matrix 
elements for the operators Tx are found in Ref. [ES 



Our coupled-channels approach is based on the Gamow state wave function for reso- 
nances. This indicates that the channel wave functions u\^^j^ni{^) have the asymptotic form of 

^ipjpnii^ipi^nii^)+'>' Fi^{kniT)) at r — i> oo for all the channels, where kni = \J2^{E — nhuj)/h^ 
is the channel wave number, and Fi and Gi are regular and irregular Coulomb wave func- 



tions, respectively 0]. This method, however, requires to solve the coupled- channels equa- 
tions in the complex energy plane and out to large distances, which is quite time consuming 
and also may be difficult to obtain accurate solutions. A much simpler alternative approach 
has been proposed, which is based on the Green's function formalism |ri] , p!^ , p!9| , p0| -|32| . This 
method was first developed for a decays by Kadmensky and his collaborators [0 and was 
recently applied to the coupled-channels problem for deformed proton emitters by Esbensen 
and Davids In this method, the coupled-channels equations are solved in the real en- 



ergy plane and the solutions are matched to the irregular Coulomb wave functions Gi^ at a 
relatively small distance Vmatch, which is outside the range of nuclear couplings. From the 
solution of the coupled-channels equations \E'^^(r, a) thus obtained, the outgoing wave func- 
tion for the resonance Gamow state is generated using the Coulomb propagator as |p!8|,p!9 



"^jmir, a) 



dr'da ( ra 



Hcoui ~l~ 



vib 



E 



17] 




coup 



y, a') 



jm\' 1 ^ ) ) 



(10) 

where r] is an infinitesimal number and Hcoui = — + Z^e^ /r is the Hamiltonian for 
the point Coulomb field. As is well known, the single particle Green's function is expressed 
in terms of the regular and the outgoing wave functions |jl8|. The asymptotic normalization 



factors A^/"'- „j then read [|T9| 



2/i 



h knI Jo 



drrFi^{knir) ( {lpjpnl)jm 



Vcoupir, a) - 



'111 



In this way, the effects of the long range Coulomb couplings outside the matching radius 
Tmatch SiYe treated perturbatively. Computing the asymptotic outgoing flux, the total decay 
width is found to be 



E 

IpjpTlI 



fi 



I 

pjpTlI I 



1 



\ jml jmi 



(12) 
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where the normahzation factor is calculated inside the outer turning point. The 

decay half-life is then defined as 

T^/2 = ^^ln2, (13) 

where Sj is the spectroscopic factor for the resonance state. If one assumes that the ground 
state of an odd-Z nucleus is a one-quasiparticle state, the spectroscopic factor Sj is iden- 
tical to the unoccupation probability for this state and is given by 5*^ = u'j in the BCS 
approximation ||TT| . 

Let us now numerically solve the coupled-channels equations for the resonant l/iii/2 
and 3si/2 states in ^^^Re as well as the 2o?3/2 state in ^^°Re. We use the real part of the 
Becchetti-Greenless optical model potential for the proton-daughter nucleus potential 



The potential depth was adjusted so as to reproduce the experimental proton decay Q 
value for each value of the dynamical deformation parameter f3 and the excitation energy 
huj of the vibrational phonon excitations in the daughter nucleus. Following Ref. 



we 



assume that the depth of the spin-orbit potential is related to that of the central potential by 
Vso = — 0.2Vo. The charge radius Rc is assumed to be the same as Rq in the nuclear potential. 



The spectroscopic factor Sj in Eq. ( p!3D is taken from Ref. []1T|. This was evaluated in the 
BCS approximation to a monopole pairing Hamiltonian for single-particle levels obtained 
with a spherical Woods-Saxon potential. 

We first discuss the effects of the quadrupole vibrational excitations. Figure 1 shows the 
dependence of the decay half-life for proton emissions from the three resonance states on 
the dynamical deformation parameter of the quadrupole mode. The excitation energy of 
the quadrupole phonon is set to be 0.8 MeV. The experimental values for the decay half-life 
are taken from Refs. and denoted by the dashed lines. The arrows are the half-lives in 
the absence of the vibrational mode. For the decays from the 3si/2 and I/111/2 states, if one 
takes into account uncertainty in the experimental Q value of the proton emission, these 
calculations for the decay half-life in the no coupling limit are within the experimental error 



bars [^. In contrast, the calculation for the 2^3/2 state is still off from the experimental data 



by about 30% even when uncertainty of the decay Q value is taken into consideration |]TT 
The results of the coupled-channels calculations are shown by the solid lines in the figure. 
One notices that the channel coupling effects significantly enhance the decay half-life for the 
2d3/2 state, while the effects are more marginal for the 3si/2 and l/iii/2 states (notice the 
difference of the scale of the vertical axis in the figure). If one assumes that the vibrational 
properties are identical between the core nuclei ^^^W and ^^°W, the dynamical deformation 
parameter of the range 0.18 < P2 ^ 0.23 simultaneously reproduces the measured decay 
half-life for the 2^3/2, 3si/2 and l/iii/2 states. 

Figure 2 shows the results for different value of the excitation energy of the vibrational 
phonon. For each panel, with the decreasing order, the four curves denote the decay half-life 
for huj=0.6, 0.7,0.8, and 0.9 MeV, respectively. The half-life for the 3si/2 state is relatively 
sensitive to the phonon energy. For the phonon energy smaller than about 0.6 MeV, the 
curve raises too quickly as a function of the dynamical deformation parameter, and the 
consistent description among the decay rates for the three states is not possible. The range 
of dynamical deformation parameter which simultaneously reproduces the measured decay 
half-life for the three resonance states is plotted in fig. 3 as a function of the excitation 
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energy h(jj2- We see that the simultaneous description is possible only when fiuj2 is larger 
than 0.59 MeV. The value of /?2 larger than 0.25 would not be acceptable for spherical nuclei, 
but the minimum value of /?2 is always below this limit for hu > 0.6 MeV. 

Lastly, we would like to discuss the dependence on the multipolarity of the phonon mode. 
Figure 4 shows the effects of octupole phonon excitations on the decay half-life from the 2^3/2 
state of ^^°Re. As an illustrative example, we take huj^=0.2 MeV, but results are qualitatively 
the same for different values of hu^. For the octupole vibration, the enhancement of the 
half-life is too small to account for the observed discrepancy between the experimental decay 
half-life and the prediction of the potential model with no coupling. As was noted before , 
the proton decay is very sensitive to the angular momentum of the proton state. The same 
seems to be true for the multipolarity of the collective excitation of the daughter nucleus. 

In summary, we have solved the coupled- channels equations to take into account the 
effects of the vibrational excitations of the daughter nucleus during the proton emission. 
Applying the formalism to the the resonant l/iii/2 and 3si/2 states in ^^^Re and the 2^3/2 
state in ^^^Re, we have found that the experimental data for the decay half-lives for these 
three states can be reproduced if the quadrupole phonon excitation with huj2 > 0.6 MeV 
and /?2 ~ 0.18 is considered. This removes the discrepancy observed before between the 
experimental data and the prediction of the optical potential model calculation for the 
decaying 2^3/2 state in this mass region without destroying the agreement for the l/iii/2 and 
3si/2 states. A similar calculation with the octupole phonon was not satisfactory. In this 
paper, we estimated the spectroscopic factor in the BCS approximation for the monopole 
pairing Hamiltonian. It would be an interesting future work to consistently compute both 
the decay rate with the coupled-channels framework and the spectroscopic factor based on 
the phonon induced pairing mechanism proposed in Ref. . 
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FIG. 1. The decay half-hves for proton emission from the 2^3/2 state of ^^'^Re and the 3si/2 and 
I/111/2 states of ^^^Re as a function of the dynamical deformation parameter /32 of the quadrupole 
vibrational excitation of the daughter nuclei. The excitation energy of the quadrupole phonon is 
set to be 0.8 MeV. The experimental data are taken from Refs. [3,4] and denoted by the dashed 
lines. The arrows indicate the results in the no coupling limit. 
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FIG. 2. The dependence of the half-hves on the value of the excitation energy of the quadrupole 
phonon. In the decreasing order, the curves are for T1UJ2 = 0.6, 0.7, 0.8, and 0.9 MeV, respectively. 
The meaning of each line is the same as in Fig. 1. 
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FIG. 3. The range of the dynamical deformation parameter f32 and the excitation energy hco 
which simultaneously reproduces the experimental decay half-life for the 2^3/2 state of ^^'^Re and 
the Ssi/2 and l/iii/2 states of -"^^^Re. 
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FIG. 4. The decay half- life for proton emission from the 2^3/2 state of ^^"^Re as a function of 
the dynamical deformation parameter of the octupole phonon excitation of the core nucleus. The 
excitation energy is set to be 0.2 MeV. 
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